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Abstract 

a : 

At lowest order in perturbation theory, the scattering matrix element for Higgs boson 
production in association with dijets displays a strong correlation in the azimuthal angle 
between the dijets, induced by the CP-properties of the Higgs Boson coupling. However, the 
phase space cuts necessary for a clean extraction of the CP-properties simultaneously induce 
large corrections from emissions of hard radiation and thus formation of extra jets. The 
' current study concerns the generalization of CP-studies using the azimuthal angle between 

dijets beyond tree-level and to events with more than just two jets. By analyzing the High 
Energy Limit of hard scattering matrix elements we arrive at a set of cuts optimized to 
enhance the correlation, while maintaining a large cross section, and an observable, which 
is very stable against higher order corrections. We contrast the description of Higgs boson 
production in association with jets at different levels: for tree-level hjj and hjjj matrix 
t-H ■ elements, for hjj matrix elements plus parton shower, and in a recent all-order framework, 

which converges to the full, all-order perturbative result in the limit of large invariant mass 
between all produced particles. 

oo : 

1 Introduction 

O ' 

One of the prime goals of experiments at the CERN Large Hadron Collider (LHC) is the search 
for the Higgs boson(s) which, within the Standard Model (SM) and many of its extensions, 
. provide direct access to the dynamics of electroweak symmetry breaking. Once discovered, the 

focus of Higgs physics will turn to the study of Higgs boson properties, like its mass, spin, CP 
parity and the strength and structure of Higgs boson couplings to heavy fermions and gauge 
bosons. 

Among the various Higgs channels at the LHC, the production of a Higgs boson in association 
with two energetic jets has emerged as particularly promising in providing information on the 
dynamics of the Higgs sector. For a SM-like Higgs boson, the vector boson fusion process, 
qq — > qqh, is expected to provide Higgs signals in the h — > W + W~ [T], h — > t + t~ [SJE], and/or 
h — > 77 m decay channels [5j|6], depending on the Higgs boson mass, and produce crucial 
information for extracting the size of Higgs boson couplings [JHS]. A second important source of 
Higgs plus 2-jet events are gluon fusion processes, such as qq — > qqh, gq — > gqh or gg — > ggh |10| . 
For favorable values of the Higgs boson mass, gluon fusion induced hjj events should be visible 
at the LHC in the h W + W~ [11 J and h t + t~ [£21113] channels. 

Even in the presence of substantial SM backgrounds, the azimuthal angle correlations of 
gluon fusion induced hjj events can be used for establishing CP properties of the interactions 
of the Higgs boson [HHT3] . In a SM-like situation, the Higgs interaction with gluons is mostly 
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mediated by a top-quark loop. For a Higgs boson which is lighter than the top quark, the 
resulting hjj cross section can be determined to good approximation by an effective Lagrangian 
of energy dimension five, which is given by |10pi4j 

where G® u denotes the gluon field strength and G a ^ u = l/2G a pa e liV(K7 its dual. yf M = m t /v is 
the SM Yukawa coupling of the Higgs boson to top quarks. The two terms result from a tth 
and a ti^tA coupling of the (pseudo)scalar Higgs, respectively, and they lead to distinctively 
different distributions of the azimuthal angle between the two jets: the CP-even hgg coupling 
produces a minimum for (ftjj = ±ir/2 while a CP-odd Agg coupling leads to minima at </>,-,■ = 
and ±7r. The azimuthal angle modulations get particularly pronounced when the two jets are 
widely separated in rapidity. Equivalent effects are expected in vector boson fusion and have 
been discussed in |15U16| for the idealized situation of parton level events at leading order (LO). 

For a more realistic simulation of hjj events, effects from multiple parton emission must be 
included, i.e. a full parton shower analysis with subsequent hadronization and jet reconstruction 
should be performed. There are two effects which may substantially reduce the azimuthal angle 
correlations predicted with LO matrix elements: (i) the emission of additional gluons between 
the widely separated jets may lead to a decorrelation of azimuthal angle distributions and (ii) 
even the very definition of the two tagging jets becomes ambiguous in the multi-jet environment 
which is expected after parton shower. 

A first analysis of hjj events with full parton shower simulation was performed by Oda- 
giri [T7j, which, however, was built on LO gg —> h hard matrix elements which do not provide 
the full dynamics of non-trivial azimuthal angle correlations of the produced partons. A subse- 
quent Herwig analysis with parton shower corrections to hard matrix elements for LO Higgs 
+ 2 and Higgs + 3 parton production [18] indicated that parton shower effects do lead to some 
deterioration of the azimuthal angle correlations of the two leading jets (defined as the highest 
transverse momentum jets). The observed deterioration was larger, however, than what was ob- 
served in a full NLO calculation of hjj production in the effective Lagrangian approximation |19] . 
This discrepancy suggests that a parton shower approach generically leads to an underestimate 
of jet azimuthal angle correlations. The most likely cause is that the parton shower generates 
radiation which is basically flat in azimuthal angle and, hence, all events, where one of the tag- 
ging jets arises purely from the parton shower, have lost the azimuthal correlation. In addition, 
hard radiation can lead to substantial changes in the angle of the selected tagging jet which 
results in further decorrelation. 

The full complexity of the jets in terms of the particles arising from the shower initiated by 
a few hard partons can only be described within a general purpose Monte Carlo [20H22j with 
a hadronization model. However, the description reached in such models by starting from soft- 
and collinear approximations may not be satisfactory for the hard (in transverse momentum), 
perturbative corrections; the description of the hard radiation can be repaired by matching 
[23H25] to fixed order calculations; but this means that the effects of such hard emissions are 
not resummed, and furthermore, virtual corrections (and the resulting weighting of samples 
with varying jet multiplicity) are estimated using only the Sudakov factor from the shower. 
When the interest is in the number of jets and their topology, rather than the description of 
the constituents of each jet, then the hard, multi-parton matrix elements can be approximated 
to all orders, based on the factorization properties of the scattering matrix element in the limit 
of large invariant mass between each hard particle, known as the High Energy Limit or Multi- 
Reg ge- Kinematics [26H28] . These ideas were developed further in Ref. [29--32J, where a partonic 
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Monte Carlo program was constructed, which captures to all orders the leading logarithmic 
behavior for emissions under large invariant mass. 

The goal of the present paper is two-fold. First, we present an improved observable, replacing 
the azimuthal angle between the leading jets as a probe for the CP-properties of the Higgs 
coupling to heavy quarks. The improvements are then assessed by simulating Higgs events in 
association with two or more jets at different levels of complexity. In Section [2] we first discuss 
how insight into the structure of the perturbative corrections in the High Energy or Multi- 
Reg ge Kinematic (MRK) limit can assist in designing observables, which extract the relevant 
kinematic information of the Higgs Boson vertex and are stable against higher order corrections. 
The resulting redefinition of the azimuthal angle between the two final state partons, in terms 
of the angle between jet clusters, is then first probed in Section [31 where we compare results 
for full fixed order matrix elements: the tree-level calculations for hjj and hjjj production are 
juxtaposed and the idealized final states in the two- and three-parton configurations are used 
for deriving upper bounds on the azimuthal angle correlations which can be expected in a full 
simulation. For this full simulation we compare two approaches in Section H) the conventional 
parton shower is generated with Herwig++ [20 while the Multi-Regge Kinematics is simulated 
with the programs developed in Refs. [29H32j . A summary and final conclusions are drawn in 
Section [5j 

2 Lessons From the High Energy Limit 

In this section we will discuss how the insight into the structure of the perturbative corrections 
in the MRK limit can assist in designing observables, which extract the relevant kinematic 
information of the Higgs Boson vertex and are stable against higher order corrections. 

The MRK limit of the 2 — > n scattering process is characterized by a large invariant mass 
between each of the produced particles, each of a fixed, perturbative transverse momentum. 
Specifically, for the scattering resulting in jets with momenta pi, • • • ,p n , the MRK limit is given 
by 

Vi G {2, . . . , n - 1} : » y { > y i+1 
V«, j ■■ \pi±\ ~ \Pj±\, \li±\ ~ hj±\. 

Here qjj_ is the transverse part of the momentum of the i th i-channel propagator, defined as 

qi = Pa-J2 l j=iPj- 

In this extreme kinematic limit, the hard scattering matrix element simplifies for two reasons: 

1. The contribution to jet production is dominated by the partonic channels which allow a 
color octet exchange between all pairs of neighboring particles in rapidity. Furthermore, 
in the MRK limit all such partonic channels have a universal behavior, dictated by the 
residues in i-channel momenta (see Ref. |3iy32j for more details). 

2. The kinematic invariants relevant for the description of the dominant part of the scattering 
amplitude simplify in the MRK limit, resulting in a dependence on transverse momenta 
only. 

Recent efforts |29H32j have concentrated on developing a formalism exploiting the universality 
of the dominance of i-channel poles (point 1), without introducing unnecessary kinematic ap- 
proximations, thus obtaining a formalism which can be applied with good results away from 
the strict MRK limit. However, the discussion in the remainder of this section will be based on 
applying the full kinematic approximations of this limit. 
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In the MRK limit, the square of the tree- level hard scattering matrix element takes the 
following form for Higgs boson production in association with n gluon jets, when the rapidity 
of the Higgs boson is in-between jets j and j + 1:- 




where q a ± = — Ya=i Pi±i where j is the number of gluons with rapidity smaller than that of the 
Higgs boson, and = q a ± — P/u_- In this limit, the contribution from quark-initiated processes 
is found by just a change of one color factor Ca — > Cf for each incoming gluon replaced by a 
quark. The effective vertex for the coupling of the Higgs boson to two off-shell gluons through 
a top loop is in the combined large-m^ and MRK limit 



c^(q a ±,qb±) = i ^ (|ph±| 2 - ha±\ 2 - hb±\ 2 ) , 

yt a s ( 4 ) 

A = ~SM o ' V = 246 GeV - 

yf lvl 67TV 

Since p h = q a - q b , C H is given by 

C H (q a ±,qfe±) = -i A |q oJ _| | cos(^ qa±qb± ) = -i A q a± ■ q b _L (5) 

In the case of a CP-odd coupling of a pseudo-scalar to the fermions in the loop, the high-energy 
factorization properties and the formula in Eq. ([3]) still hold, with C H replaced by 

C A (q a _L,<ib±) = i B |q a ±| \q.bJ-\ sin(0 qa±qi)± ) = % B (q a _L x q 6± ) • z, (6) 

where z is a unit vector along the beam axis whose sign is given by the z component of q a — qb- 
In the notation of Ref . [11] , 

yfM " 2 TTV' K ] 

The CP-properties of the coupling of the (pseudo-) scalar are reflected in the dependence on the 
azimuthal angle between q a and qi, of the effective vertex for the coupling to gluons via a fermion 
loop (e.g. C H or C A ). It is this dependence that the observables and cuts should emphasize. In 
the case of Higgs boson production in association with just two jets, q a j_, qbj_ are (up to a sign) 
given by the transverse momentum of the jets (at lowest order in perturbation theory), and 
thus the CP-properties of the Higgs boson coupling can be extracted by studying the azimuthal 
angle between the two jets [IB] . 

Eq. ([3]) is valid for the rapidities of the Higgs boson within the rapidity interval spanned 
by the partons. In fact, the factorization of the amplitude into a part which depends on the 
momenta of the partons with rapidities much smaller than that of the Higgs boson, a Higgs boson 
vertex, and a function depending on the momenta of partons with rapidities much larger than 
that of the Higgs boson is valid also when no requirement is placed on the invariant mass between 
the partons on either side of the Higgs boson [34] . Therefore, we will divide the observable jets 
into two groups and require that all jets are well separated from the Higgs boson direction, 

3ja : Vj a < Vh, 3j b : y jb > y h , 
Vj e {jets} : \yj - y h \ > y sep 
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In the following we always use the kt jet algorithm as implemented in Ref. [35] with R = 0.6 
and investigate jets with a transverse momentum greater than 40 GeV. 

While it is clear that in LO hjj production the CP-structure of the Higgs boson vertex 
can be revealed by studying the azimuthal angle between the two jets, the question arises which 
angle to study in events with more than two jets. Traditionally the observable applied to extract 
the CP structure has been the azimuthal angle between the hardest jets of the event pT t fTB l fTH]. 
within a set of cuts based on the rapidity separation between the two hardest jets. However, 
Eqs. ©-([5]) suggest that the CP structure is more clearly revealed if instead the azimuthal angle 
between q a and (ft is studied, with 

la = Yl Pi> qb = Yl p i> 

je{jets :2/j <y h } je{)ets: yj >y h } (9) 

^2 = ^(q a ±,qb±) 

In the MRK limit, (1/cr da/dfa) is stable against higher order corrections, in so far as soft 
radiation, which is not picked up as jets, does not impact the angular distribution. 



3 Fixed Order Studies 

3.1 Tree level studies for hjj production 

We start by studying the tree-level predictions for hjj production for a 14 TeV proton-proton 
collider and a Higgs boson mass of 120 GeV. We use the central pdf set from MSTW2008 [55] . 
and initially we use a minimum set of cuts: 

p i± >40GeV, y ja <yh<yjb, \yj aijb \ < 4.5. (10) 

The main conclusions of this study are insensitive to the choice of renormalization and factoriza- 
tion scale. In the /ijj-analysis we set aj. — > a s (pj a ±)a s (pjb±)a1(mH), with factorization scales 
for the pdfs of fj,f ia = Pj a ±, /-*/,& = Pjb±- The matrix elements are extracted from MadGraph [37 , 
where the mj — > oo limit is used. 

The distribution of the cross section wrt. azimuthal angle between the two jets is often 
distilled into a single number 

A = v(\0jah\ < tt/4) - q(7r/4 < \</> jajb \ < 3tt/4) + a(\(f> ja j b \ > 3tt/4) 
°{\<t>j a h\ < V4) + d(vr/4 < \ct> jajb \ < 3tt/4) + <r(\<t> jajb \ > 3vr/4) ' 1 ' 

This obviously gives — 1 < Aa < 1, with A^ = representing no azimuthal correlation between 
the jets. If there are no other sources of angular dependence, then a CP-even coupling of the 
Higgs boson to two off-shell gluons through a top-triangle leads to a positive value for A^, 
whereas a CP-odd coupling results in a negative value for A^. 

In Fig. [T] we plot the value for A^ as a function of Ay and y* given by 

Ay=\y ja -y Jb \, y* = ^ - «k±Ik. (12) 

The bin-size is .25 x .25 units of rapidity. The insert shows the cross section per bin. Because of 
the requirement of the Higgs boson to be in-between the jets in rapidity, we have < y* < Ay/2. 
We see that A^ is largest when the jets are separated by a large rapidity interval (Ay large), and 
the Higgs boson is produced near the rapidity-center of the two jets (y* small). This is caused 
by two effects: 
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Figure 1: vs. Ay and y* with a bin size of .25 x .25 units of rapidity, for a SM Higgs boson 
of mass m/i = 120 GeV. Insert: da j (dy* dAy) in fb/bin. 

1. This kinematic setup approaches the MRK limit, where the matrix element for the gluon 
channels have the same kinematic dependence as the quark-initiated ones. This depen- 
dence is induced by the singularity from a gluon exchanged in the t-channel, coupling to 
the Higgs boson field through a top triangle. The effects of other gluon-Higgs couplings 
through top loops are suppressed, and only the one giving rise to the strong azimuthal 
correlation survives (see Ref. j29ft32j for further details). 

2. The azimuthal angular dependence of the t-channel propagators is suppressed, and the 
only azimuthal dependence left is induced by the tensor structure of <7*<7*/i-coupling which 
leads to Eq. ©. 

The insert in Fig. [1] shows that the cross section is peaked at Ay ~ 1 and decreases for increasing 
Ay and y*. On Fig. Q]we have also indicated the limits in the (Ay, y*)-plane corresponding to 
the three values y sep = 0.5,0.75, 1, i.e. the lines correspond to 

mm(\y h - y ja \,\y h - Vj b \) = y sep • (13) 

Good analyzing power (large A^) and a large rate are simultaneously retained by imposing a 
y S ep cut only (as opposed to Ay and/or y* cuts). Thus we require that the Higgs boson be 
produced between the jets, and that the minimum rapidity distance between the Higgs boson 
and the two jets be larger than a minimal value, y S ep- 
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Figure 2: (1/cr da/dfyj) for /yj-production at tree- level with the cuts of Eq. (|14p and y sep = 
0.5,0.75,1.0. 

Based on the findings presented in Fig. [JJ and the structure of the re-jet scattering amplitude 
in the MRK limit given by Eq. ([3]), we suggest applying the following cut in Higgs boson+multi- 
jet samples, when the emphasis is on extracting the CP-structure of the Higgs boson coupling. 

p j± > 40GeV, 3j a Jb G {jets} : y ja < Vh < Vjb, Vj G {jets} : \y h - yj\ > y sep , \yj\ < 4.5. (14) 

In Fig. [2] we plot (l/ada/dcpjj), where (fijj is the azimuthal angle between the two jets (defined 
with the sign according to Ref. [16J), with the cuts of Eq. (fT4"j) and for y scp = 0.5,0.75,1.0. 
We see that as y sep is increased, the increase in noted in Fig. [1] is reflected in the angular 
distribution approaching that of cos 2 (<?!>), as expected from Eq. ([3{lj. The cross section and values 
of Aff, for the series of cuts are summarized in For comparison we mention that with the same 
scale choices, but using the standard weak boson fusion (WBF) cuts, the LO contribution to hjj 
through gluon fusion is 230 fb with A^ = 0.46 (see e.g. Ref. [30] for further details). Table [TJ So 



Usep 


a hjj 


A^ 


0.5 


743 fb 


0.36 


0.75 


553 fb 


0.41 


1.0 


403 fb 


0.45 



Table 1: The value for the tree- level hj j-cross section with the cuts of Eq. (|14|) and three values 
of y scp - 

1 Obviously, the cos 2 (^>)-behavior of the cross section is reached only once not just the numerator is well 
approximated by that of Eq. §3$, but also the £-channel propagator momenta have lost dependence on all but the 
transverse degrees of freedom. This is true in the MRK limit and assumed in Eq. ([3J, but the sub-asymptotic 
effects are incorporated in the formalism developed in Ref. |29H32] 
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Figure 3: (1/<t dcr/d<f)jj) for production of the CP-odd higgs boson A plus two jets at tree-level 
with the cuts of Eq. (fT4"|) and y sep = 0.5, 0.75, 1.0. The scales were chosen as described in Section 
IP 

the cuts in Eq. (|14| with y sep = 1 display the same strong azimuthal correlation while almost 
doubling the cross section. 

In Table [2] we list the cross sections and values of found for the case of a CP-odd 
coupling of SM-strength to a pseudo-scalar through the effective vertex of Eq. ([6]). We note 
that the magnitude of A^ is almost unchanged compared to the CP-even case, just the sign 
is reversed. In the remainder of this study we will focus on the case of the CP-even Standard 



2/sep 




A^ 


0.5 


1730 fb 


-0.36 


0.75 


1271 fb 


-0.41 


1.0 


908 fb 


-0.46 



Table 2: The value for the tree-level Ajj-cioss section with the cuts of Eq. (|14|) and three values 
of y sep . The scales were chosen as in Section [4JJ 

Model Higgs boson. However, it is clear that the cuts and analysis developed in the following 
will be well suited for the extraction of the CP-parameters, even in the case of a mixture of 
CP-even and CP-odd couplings. 

3.2 Extracting the CP-Structure with a Fit 

In experimental data, there will of course be a non-negligible background to the gluon fusion 
/y'j-signal. While the background is predominantly flat in <j> [11H13] . a flat pedestal in the 
distribution (da /deft) does not leave invariant. We will therefore introduce a different measure 
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of the azimuthal correlation by fitting the distribution to the following form 

f(4>) = A + |p cos^ + D) + C<P\ (15) 

where the parameters A, B, C, D are fitted. The </> 4 -term is found to describe well the impact of 
the sub-asymptotic angular dependence. In the true asymptotic form of the scattering with a 
CP-even coupling and no background contribution, A = %B, D = C = 0. For C = D = one 
has A^ = B I A. A CP-odd coupling would set D = tt/2. The size of B (measured in fb/rad) 
is a direct measure of the ability of the cuts and definition of angle in resolving the azimuthal 
correlation induced by the CP-nature of the Higgs boson coupling. Since in the current study 
we are interested in the SM Higgs couplings only, all the fits return a value for D consistent with 
0. The quantity B/A measures the quality of the given sample for extracting the CP-properties. 
The values for B/A obtained in the fits to the hjj tree- level distributions for y se p = 0.5,0.75,1 
are 0.32, 0.38, 0.44, respectively. We notice that as expected, the agreement between A^ of 
Tableland B/A improves for increasing y sep . 

3.3 Higgs Boson Production in Association with Three Jets 

We will now discuss the stability of the shape of (1/cr da/d(f>) found for hjj at tree level against 
real emission higher order corrections within the cuts of Eq. (|14|) . when <f> is defined according 
to Eq. ([9]). The results of the MRK limit indicate that (1/cr da/d<f)) should be stable, when is 
defined according to Eq. ([§]), rather than just the azimuthal angle between the two hardest jets. 
The results for the tree-level hjj j-cross section and the extracted A^ are listed in Table [3l We 
note that with the definition of the azimuthal angle of the multi-jet events according to Eq. ([9]), 
for each value of y sep , the extracted value for A^ changes only very little compared to the values 
extracted from the two-jet case. For comparison, we note that if the WBF cuts are applied to 
the two hardest jet in the three-jet sample, the tree-level cross section for hjjj is found to be 
just 76 fb, with A^ = 0.37 (a change of 0.09 units), when <f> is the azimuthal angle between the 
two hardest jets. 

The differences in the distribution (1/cr da/d(j)) (with c/> defined according to Eq. ©) in going 
from the two to three-jet production is less than the differences between each neighboring curve 
in Fig. [21 corresponding to varying y sep . The stability is a direct result of the definition of the 
angle according to Eq. Q , designed from the insight obtained by analyzing the MRK limit of the 
amplitudes, and optimized to reflect directly the CP-structure of the Higgs boson coupling. This 
stability in the reconstruction and extraction of the azimuthal angular dependence is encouraging 
for the possibility of determining the CP-properties of the Higgs boson coupling from LHC data. 

In the Section 0] we will inspect the stability of this azimuthal angular dependence in two 
models for all-order perturbative corrections: 1) the Parton Shower model as implemented in 
HerwigH — h, supplemented with the tree-level /ijj-amplitudes as described in Ref. [18J, 2) a 
model for both the real and virtual corrections, based on the MRK limit of amplitudes, as 



?/sep 


a hjjj 


A^ 


B/A 


0.5 


365 fb 


0.34 


0.30 


0.75 


232 fb 


0.39 


0.34 


1.0 


150 fb 


0.42 


0.41 



Table 3: Values for the tree-level hjj j-cross section, A^ and the fitted B/A ratios with the cuts 
of Eq. (fl~4"j) and three values of y sep (4> defined according to Eq. ((9|)). 
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discussed in Refs. [29l[30]. These studies will probe the sensitivity to the azimuthal correlation 
against soft radiation, which is not picked up in the visible jets, and (especially in the second 
case) further hard radiation. We will choose y sep = 0.75 - the specific value of y sep does not affect 
the conclusions beyond the variation indicated in Tables [T][3J Larger values of the accepted cross 
section could be obtained by choosing a smaller value for y sep . 

In the remainder of this section we will test two further predictions on the azimuthal distri- 
bution obtained from limiting behavior in Eq. ([3|) of the full /ijjj-amplitude. According to the 
cuts of Eq. (|8|) , the three-jet events will have two jets on one side (in rapidity) of the Higgs boson, 
and one jet on the other. We will call the two jets on the same side j a ,jb, ordered according to 
hardness. The single jet on the other side of the Higgs boson is called j c . First (in Section [3.3.ip . 
we see that the distribution in the azimuthal angle defined according to Eq. for three-jet 
events should be independent of e.g. the angle between j a ,jb, as long as the azimuthal angle 
of the sum of these two jet vectors is fixed. Secondly, (in Section I3.3.2|) . we investigate the 
prediction that the strong azimuthal correlation is displayed also in events where the jets are 
ordered in transverse momentum as follows: p±j a > P±,j b > P±,j c > 40 GeV, i.e. the two hardest 
jets are on the same side in rapidity of the Higgs boson, and the softest jet is on the other side. 

3.3.1 Dependence on the Azimuthal Orientation of Two Jets 

In Fig. H] we have plotted the tree- level predictions for (1/ a da /d(p) for a sample of Higgs boson 
plus three-jet events with-in the cuts of Eq. (y sep = .75), for three bins of the azimuthal 
angle between the two jets on one side of the Higgs boson: 

51 ■ 4>i a j b < vr/4 

52 : vr/4 < <f> jajb < 3vr/4 (16) 

53 : 37T/4 < ioA < 7T. 

For reference, we have also plotted the distribution in the azimuthal angle in the two-jet sample. 
It is clear that the azimuthal dependence with the definition of Eq. Q for the three-jet samples 
is extremely stable when going to the three-jet samples SI and S2. Only for the sample S3 is 
there a slight change - this is because this sample is forced to have the two jets on one side of the 
Higgs boson back-to-back. Since the cross section is dominated by configurations where the jet 
transverse momentum is close to the lower cut-off, the sum of the transverse momenta of the two 
jets will typically be close to 0, and the square of the momentum of the t-channel propagator is 
not dominated by its transverse momentum. In other words, the last requirement of the MRK 
limit in Eq. ([2]) is violated in this sample, introducing a large dependence of the full t-channel 
propagator on the azimuthal orientation, so the azimuthal dependence is not described by the 
numerator alone. This is one reason why the formalism of Ref. |29H32| (which implements 
the full momentum-dependence on the t-channel propagators) offers an improvement on the 
BFKL formalism |26 t l28 l l3"5]. Fig. H] also shows the fit according to Eq. (|15|) to each of the 
four distributions. It is clear that the functional form in Eq. (|15p describes the distributions 
very well. The value for B/A obtained in the fit are 0.39,0.35,0.24 for the samples <S1, .S2, .S3, 
respectively, compared with 0.34 for the full three-jet sample with the same y sep = 0.75. 

3.3.2 Dependence on the Hardness of the Three Jets 

Let us now turn to the case with two harder jets on one side in rapidity of the Higgs boson, 
and a softer on the other side, all above 40 GeV in transverse momentum. The reason for 
studying these event configurations is to determine (using the full tree-level hard scattering 
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Figure 4: The angular distribution (according to Eq. Q) for y sep = 0.75 for both the two-jet 
sample and three different sub-samples of the three-jet sample (all jets above 40 GeV in trans- 
verse momentum — see text for further details). The cuts and the definition of azimuthal angle 
in samples of more than two jets ensure that the distribution is very stable against perturbative 
corrections. The lines are fits on the form of Eq. (|15p to the histograms. 



matrix element) whether in this skewed setup the scattering displays the angular correlation 
expected from the analysis in Section [2j Specifically, we choose the following cuts 

P±,ja > P±,j b > P±,jc > 40 GeV > 

Vjajb <Vh< Vj c or y jaJb > y h > y jc (17) 
\Vh -Vj\ > ysc P , \yj\ < 4.5. 

With the standard scale-choices, the cross section passing these cuts (with y sep = 0.75) is 82 fb, 
i.e. about 1/3 of the three-jet rate without the additional requirement of a special configuration 
of the hard jets. Note that this is exactly the fraction expected from the MRK analysis of Sec. [2j 
There is no preference for any of the three jets in the event to be the softest, so in 1/3 of the 
events in the inclusive 3-jet sample it is the jet which is separated in rapidity from the two others 
by the Higgs boson. 

In Fig. [5]we have plotted (1/cr da/defi) for hjj, and two different azimuthal distributions for 
the sample of three-jet events with these selection cuts. It is immediately clear that 1) the two 
hardest jets do not display even a hint of the azimuthal correlation of the two-jet sample (the 
only slight enhancement is when the two jets are pointing in the same azimuthal direction, with 
the collinear divergence regularized by the jet measure {kt — jet, R = 0.6), 2) the vector sum of 
Eq. still displays the azimuthal correlation induced by the Higgs boson vertex, with a value 
of B/A from the fit to Eq. (|15f> of 0.46. Therefore, if one is applying the observable of Eq. (|9|) 
in the analysis, then the events with the two hardest jets on one side of the Higgs boson can be 
included and contribute to the extraction of the Higgs boson properties. This idea is backed by 
the insight coming from the high energy factorization of the hard scattering matrix element. 
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Figure 5: The angular distribution (according to Eq. Q) for y sep = 0.75 for both the two-jet 
sample and a three-jet (all above 40 GeV in transverse momentum) sample of the two hardest jets 
on one side of the Higgs and a softer on the other (see text for further details) . The characteristic 
azimuthal dependence is not reflected in the azimuthal angle between the two hardest jets. This 
dependence, induced by the CP-properties of the Higgs boson couplings, is displayed however, 
when the sum of jet vectors according to Eq. Q is considered. The lines are fits on the form of 
Eq. (fT5|) to the histograms. 

4 Correlations in a Multi-Jet Environment 
4.1 Parton Shower 

In order to obtain results which include higher order emissions, we used the parton level Monte 
Carlo VBFNLO [39J to calculate gluon fusion processes with hjj signature, generating Les 
Houches Event Files, which were then used as an input for the event generator Herwig++ 
[20] . The renormalization scale was set to a* = u s {pj a ±)a s {pjbi_)a 2 s {rnH) and the factorization 
scale was chosen as fif = y/Pj a ±Pjb±- The computations were done using CTEQ6L1 parton 
distribution functions |40| . which is different from the choice of Section [3l However, this has 
only minor effects and leaves the main conclusions of our findings untouched. 
We generated unweighted events using the following cuts at the parton level, 

p j± > 30 GeV, \ yj \ < 7, \Rjj\ > 0.3, (18) 

and then applied the cuts of Eq. (|14p on the recombined jets (which were formed using the kt 
jet algorithm with a cone parameter of R = 0.6 as in the previous section) after the shower, 
hadronization and decays. An underlying event was not simulated in our studies. The choice of 
rather weak cuts at the matrix element level leaves the parton shower evolution a great latitude, 
especially concerning effects altering rapidities of the jets. 

The HerwigH — h parton shower has a built in veto to forbid too hard shower emissions. The 
default value for the veto scale is the highest transverse momentum of the outgoing partons 
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Figure 6: (1/cr da/d(j)) for the hjj matrix element calculation and the results after parton shower, 
both with the cuts of Eq. (fT4"|) and y sep = 0.75. <f> is defined according to Eq. Q. The numbers 
for the azimuthal angle asymmetry measure were extracted from the fitted function. 

in each event at the matrix element level. This had to be lowered to the smallest transverse 
momentum as, otherwise, the shower produced too many hard emissions which overcame the 
characteristics of the underlying process. 

4.1.1 Higgs Boson Production in Association with Two Jets 

Fig. [6] shows the results after the parton shower simulation for y scp = 0.75. As expected, 
additional emissions created by the shower lead to a decreasing amplitude of the sinusoidal 
modulation, from B/A = 0.38 at hjj matrix element level to B/A = 0.26. However, the 
characteristics of the azimuthal angle distribution are not significantly affected by the shower. 
Events with exactly 2 jets of pt > 40 GeV constitute 73% of the cross section. In this exclusive 
2-jet sample, the shower is expected to better preserve the directions of jets which were already 
present in the matrix element calculation. Indeed, for this sub-sample, the fitted azimuthal 
asymmetry decreases by a somewhat smaller amount, to B/A = 0.30. 

Since the contribution from events with more than 2 jets is modest, it is not surprising that 
the angular distribution in the 2-jet inclusive sample closely resembles the tree level result. 

4.1.2 Higgs Boson Production in Association with Three Jets 

The situation changes, once the parton shower is asked exclusively for events with three jets in 
the final state which, however, are still based on the hard hjj matrix elements. In Fig. [7J the 
azimuthal angle distribution for these events is plotted as the red dotted histogram and compared 
to the corresponding distribution at the hjjj hard matrix element level (black dashed histogram) , 
which was already shown as the red curve in Fig. [5j The difference between the results from the 
parton shower and that of the fixed order calculation is larger than in the exclusive 2-jet case. 
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Figure 7: The azimuthal angle correlations for the hjjj sample with transverse momentum of 
each jet above 40 GeV and y sep = 0.75. The figure shows the result of the matrix element 
calculation and the distribution after the parton shower simulation based on the hjj matrix 
element calculation. 

As the shower includes effects which go beyond the tree- level matrix element calculation, some 
additional decorrelation is expected. However, these effects of the radiation beyond that of the 
tree-level are small, as shown in the previous section for the exclusive two-jet sample. Rather, 
as we will show next, a substantial part of the decorrelation is due to phase space regions where 
the parton shower completely fails to produce the azimuthal angle modulation inherent in the 
exact matrix elements. The shower algorithm currently implemented in HERWIG++ uses only 
azimuthally averaged emissions, and we will in the following check the description of the three-jet 
events, where at least one jet must arise from the parton shower description. 

Let us investigate in detail the kinematically restricted case from Section I3.3.2I with the two 
hardest jets j a and % on one side of the Higgs Boson and the third hardest jet j c on the other 
side. As noted in Section I3.3.21 this event configuration corresponds to roughly a third of the 
full three-jet sample in the case of the tree-level calculation. We focus on this sample, not 
because the parton shower emission generates a good description of the full three-jet sample, 
but because the failure in this particular sub-sample is so spectacular and makes it easy to 
draw conclusions. We demand p±j a bc > 40 GeV and the rapidity cuts of Eq. (|14j) . The results 
are shown in Fig. Ela) as the red dash-dotted histogram, while the black dashed histogram 
depicts the hjjj fixed order calculation on parton level. Unlike the cases which were examined 
before, the azimuthal characteristics of the parton level appear to be completely lost and the 
fit returns a value of B/A = 0.01. The analysis was repeated with the transverse momentum of 
the softest jet lowered: p±j a b > 40GeV, p±j c > 20GeV, in order to investigate if this hierarchy 
in transverse scales would help the description obtained within the parton shower. However, 
this did not change the azimuthal distributions noticeably, neither for the results obtained with 
tree-level hjjj matrix elements, nor for the calculation with /ij'j+parton shower. 
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Fi gure 8: The angular distribution for y se p — 0.75 for the sample with two hard jets on one side 
of the Higgs Boson and a softer (still with p± > 40GeV) on the other. The two graphs in the 
upper plot depict the results of the matrix element calculation and after parton shower. The 
lower plot shows the observable for two different subsamples: One where the events pass the cut 
of Eq. ([5D already on matrix element level with just two partons (denoted with 'PL passed') and 
one where they do not ('PL failed'). The black dotted curve is the same as in the upper plot, 
but normalized to the cross section of the subsample which passes the cuts at the /ijj-parton 
level. 
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Figure 9: Absolute values of rapidity differences for y sep = 0.75 for the sample with two hard 
jets on one side of the Higgs Boson and a softer (with pj_ > 40GeV) on the other. The blue 
dashed graph shows the rapidity difference between the two clusters of jets on each side of the 
higgs boson, after parton shower. The black dot-dashed graph shows the rapidity difference of 
the two partons at the matrix element level for the same sample of events. Again, this is split 
into two different subsamples, one where the events pass the cut of Eq. (JSj) at the parton level, 
and one where they do not. 



The question why the parton shower approximation fails to describe the hjjj matrix ele- 
ment in this special setup can be answered by having a closer look at the hard matrix element 
description underlying the events. One would expect the description from the parton shower to 
agree better with the one arrived at using the full three-jet matrix element in cases where the 
partonic /yj-configurations already fulfil the cuts in Eq. ()14[) and the two same side jets from 
the shower can be produced by final state collinear splittings. This is indeed the shown 
by the green dashed histogram in Fig. [8](b) . These events alone yield B/A = 0.23 while their 
complement, where the original hjj parton configurations fail the cuts of Eq. (|14p . have lost any 
azimuthal angle correlation and produce B/A = —0.07. 

The clearest difference between events which pass at hjj parton level and those that fail can 
be seen by comparing rapidity separations of partons and produced clusters of jets. Fig. [9] shows 
the rapidity difference 

|Ay| = |y Ca -y C6 | (19) 

of the two jet clusters c a and c^, which are present after the recombination of hadronic final 
states, as the blue dashed line. For the same sample of events, the rapidity difference of the two 
partons at matrix element level is plotted via the black dash-dotted histogram for all events, and 
this curve is split into events wich fail the cuts at parton level (red solid histogram) and those 
that pass (green dotted histogram). The sample with negligible azimuthal angle correlation in 
Fig. [8](b) was evolved from events which feature only a small rapidity separation between the 
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two outgoing partons on matrix element level and which, therefore, have both partons on the 
same side of the Higgs boson. The parton shower does not generate the azimuthal correlation 
of the clusters on opposite sides of the Higgs for such events. 

As a conclusion, the question of azimuthal angle decorrelation in the setup with two hard 
jets on one side of the Higgs Boson and a softer one on the other is too exclusive to be handled 
reliably within the parton shower approximation presently implemented in Herwig++: The 
azimuthal angle averaged parton shower is not the right tool for the observable in this case. It 
is likely, however, that the description would be improved if spin-correlations were implemented 
in both the parton shower |41U42j and the interface to the matrix element. 

4.2 Multi Parton Emissions from High Energy Factorization 

In this section we investigate the description of the production of a Higgs boson in associ- 
ation with multiple jets, as described by the formalism developed in Ref. [291130], with the 
improvements of Ref. |3iy32j implemented. This results in an all-order partonic estimate of the 
production of a Higgs boson plus at least two jets, which is exact in the limit of large invariant 
mass between all (hard) partons, i.e. specifically in the limit where the real radiation produces 
additional hard jets. The calculation is completely exclusive in the momenta of all produced 
particles. We can therefore perform a detailed investigation of e.g. the jet activity in the events. 

The all-order calculation results in a cross section within the cuts of Eq. (|14p of 456 fb for 
Usep = 0.75. This is less than the tree- level rate (listed in Table [1]). There is a reduction in the 
hjj tree-level cross section as a direct result of the veto of hard jets in the region of 0.75 units 
of rapidity either side of the Higgs boson. The exclusive jet rates found within this calculation 
are plotted in Fig. [101 

The exclusive hard two-jet rate amounts to just roughly 60% of the total rate within these 
cuts, as compared to 73% found in the study of Section 0J] based on a parton shower of partonic 
hjj-st&tes. It is therefore beneficial to ensure a solid extraction of the CP-structure of the Higgs 
coupling also in events with strictly more than two jets, as proposed in this study. In Fig. [11] 
we compare the differential distribution in the azimuthal variable defined in Eq. §§§ for the 
tree- level hj j-calculation, and the all-order hj j-calculation of Ref. |29H32| . The stability of this 
distribution against higher order corrections as implemented in Ref. |29H32j is partly a result of 
the assumption of a dominance of the f-channel pole, but this dominance is directly verified by 
comparing the results of the formalism order-by-order. However, the strong correlation seen at 
the tree-level (as a result of the dominance of the t-channel pole) could have been spoilt in the 
all-order calculation by all the partons which are not included in the momenta of the hard jets, 
on which the azimuthal observable operates. However, as we see in Fig. [TT1 operating only on 
the momenta of the jets with a transverse momenta larger than 40 GeV still allows for a very 
clear extraction of the CP-structure, when the observable of Eq. (|9|) is used. 

Finally, in Fig. [12] we compare the description of (exclusive) three-jet events between tree- 
level and the resummation of Ref. [29H32] . Specifically, we compare the description within the 
special event sample of three hard jets above 40 GeV transverse momentum, but with the two 
hardest on one side, and the softest at the other side (in rapidity) of the Higgs boson. 

While the description of this exclusive three-jet sample differed between the description 
based on the full hjjj matrix element and that based on /ijj+parton shower, we see that the 
description based on Ref. |29H32j gives results which are very similar to those obtained in the 
calculation based on the full three-jet tree-level matrix element. This is a testament to both 
the quality of the all-order approximation applied, and the stability of the angular observable, 
defined according to Eq. Q, against higher order corrections. 
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Figure 10: The exclusive jet rates within the cuts of Eq. ([8]) as obtained in the all-order calcu- 
lation of Ref. [29 32j . All jets have a transverse momentum greater than 40 GeV. 




Figure 11: A comparison between the results obtained for inclusive hjj at tree- level and for 
the all-order approximation of Ref. [29H32] for the azimuthal correlation defined according to 
Eq. ([9|). The similarity is a testament to both the quality of the approximations in Ref. [29H32] 
and the stability of the observable against higher order corrections. The lines are fits on the 
form of Eq. (j 15|) to the histograms. 
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Figure 12: A comparison between the results for the azimuthal correlation defined according to 
Eq. (J9j) obtained for the exclusive hjjj-v&te with the transverse momentum of the one jet at one 
side of the Higgs boson smaller than that of either jets at the other side of the Higgs boson, at 
tree-level and for the all-order approximation of Ref. |29H32] . The similarity is a testament to 
both the quality of the approximations in Ref. [29H32] and the stability of the observable against 
higher order corrections. The lines are fits on the form of Eq. (|15p to the histograms. 



5 Summary and Conclusions 

Based on the analysis of the Multi- Reg ge- Kinematic (MRK) limit of the tree-level matrix ele- 
ments for Higgs boson production in association with multiple jets, we optimize the extraction 
of the CP-properties of the Higgs boson coupling to heavy fermions in two steps. First, we 
have designed a set of cuts which enhances the sought-after behavior. Secondly, we generalize 
the azimuthal angle observable used to extract the CP-properties of the couplings in Higgs 
boson plus dijet events to the case of three or more jets. This is desirable, since roughly half 
of the inclusive Higgs-boson plus dijet sample will contain more than two jets (for a 14 TeV 
proton-proton collider and the cuts employed in the present analysis). 

We have thoroughly tested the predictions obtained from the MRK limit against calculations 
using full tree-level matrix elements, and find that the expectations from the MRK limit are 
all respected by the full tree-level calculations, and that the constructed observable leads to an 
extraction of the CP-properties of the Higgs-boson couplings which is extremely stable against 
perturbative corrections. 

We have compared the predictions for both Higgs boson plus two jet and for three-jet events 
obtained in two all-order perturbative approximations against the full tree-level calculation. We 
find that the calculation based on a parton shower evolution of the hard /yj-matrix element, 
while giving a good description in large areas of phase space, fails to give a satisfactory de- 
scription of the three-jet states in those phase space regions, where the direct correspondence 
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of hard partons at matrix element level with the leading jets is lost. These phase space regions 
are numerically important. In particular, we discussed the failings in the parton-shower descrip- 
tion of the azimuthal correlations observed in the full hard matrix element. Contrary to this, 
the all-order perturbative calculation based on the formalism developed in Ref. [29H32] leads 
to results which are in good agreement with both the exclusive two and three-parton samples 
checked in this study, while giving an estimate of the sensitivity of the observables to corrections 
beyond tree-level. Again, the observable constructed in Sec. [2] is shown to lead to an extraction 
of the CP-properties of the Higgs boson coupling which is extremely stable against higher order 
corrections. 

Based on the analysis of hjj and hjjj at fixed order, and the MRK-inspired all-order resum- 
mation, we conclude that, with the right choice of observable, the analyzing power of azimuthal 
jet correlations in hjj events for the extraction of Higgs boson CP-properties can be preserved 
when going from a LO analysis to a more complete and sufficiently detailed QCD description 
of additional parton emission. While the contribution from hard radiation beyond tree-level is 
significant and can jeopardize the extraction of the CP-properties, the impact of such radiation 
is minimized by the use of the observable advocated in Eq. ©, which depends on all the jet 
vectors of the event. In particular, also the NLO QCD calculation of Ref. [19J should be able 
to provide reliable predictions for the extraction of CP-properties in gluon fusion induced hjj 
events. Stronger decorrelation effects, which were observed in the past when performing parton 
shower simulations of hjj events |17pi8j are, to a considerable extent, due to present limitations 
of the parton shower programs and the use of azimuthal observables, which are unstable against 
radiative corrections. 

The current study also provides a clear demonstration of how the insight gained from the 
structure of scattering matrix elements in the MRK limit can assist the design of analyses at 
the LHC. 
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